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ABSTRACT 


rFOinting error standard deviations for two theodolit 


iD 


S, 
the Wild T-2 and Odom Aztrac, were Jietermined under condi- 
tions closely approximating those of range-azinuth or 
azimuth-azimuth hydrographic surveys. Pointing errors found 
for both instrumerts wer2 about 1.3 meters, and were inde- 
pendent of distance. No statistical difference between the 
errors of the two instruments was found. The accuracy of 
the interpolation methods used by the National Ocean Service 
(NOS) for range-azimuth positioning were investiyated, and 
an average inverse distance of about 2.5 meters was observed 
between interpolated positions and corresponding observed 
positions. The overall range-azimuth position errors of the 
two theodolites were then compared to positioning standards 
of NOS and the International Hydrographic Organizaticn, 
using assumed ranging standard deviations of 1.0 and 3.9 
meters. Both instruments met all standards except the NOS 
range-azimuth standard for 1:5,000 scale surveys. 
Interpolated positions may fail to meet more of the stan- 


2 
dards because of additional inherent error. 





Lage 


12 er 


IV. 


TABLE OF CONTENTS 


PSEC OMetINCIN = s 6 6 © « «@ © 6 6 « « 


A. TRE RANGE-AZIMUTH POSITIINING 


METHOD . 


Be. HYDROGRAPHIC POSITION ERROR STANDARDS 
er OD enn VSS: 5 © 6 6 «© 6 « 6 © @. 


PeeOreLNDIZES AND RANGE-AZIMUTH GEOMErRY . 
we DET ENTREGHS & e e e e e eo 2 e e 


1. 
on 
Be 


EOIN GesS. «6 4 Sw eS ow Ok 
Siew oaetCrrmrhOLS 4 «6 % sis 


Random Errors ees er isl der Ses 


cee SAM OUMEYSLONAL ERROR FISURES  . 


im 
2. 
Se 
4. 


Senechis Cc Ara ESecsatric G=ecn 


feicgemanOr Ellipse 6 ws vole 1% 
Root Mean Square Distance . 


Seazeuile> Standard Excar.-o.. 


eee ee ern On OF AN INTSRPOLATED FIX 


Uc 
2. 


ites mero. . On ALgOorz eh ms ., 


BeEOE MET OUSO St 26M, a « 2 b-< os 


EXPESIN@NT DESIGN AND IMPLEYENTATIORN 
Pepe oore OO WORK 6 9s = 6 6 «© «ws © «© « 
Be FHEODOLITE POINTING ERROR «..=. 
C. INIFERPOLATION ALGORITHM EVALUATION 


Demeenorc Ss OF EXPERIMENTAL CONDITIONS 


ReoU LTS 


AND DATA ANALYSIS .... . 


Ape Daten PROCESSING SYSTEM . 5 «2 « « 
B. POINTING ERROR DETERMINATION. . 
C. INI ERPOLATION EVALUATION ... . 


DreeoimeystS OF FACTORS AFFECTING THE RESULTS 


() 


(t 


\O 


ie 
14 


20 
20 
Ze 
AS 
24 
26 
Ze 
23 
SHS) 
Eye 
Bul 
37 
58 


43 
43 
47 
49 
51 


56 
5116) 
56 
64 
66 








Peet oe tent TON TO POSTTEOGN ERROR STANDARDS . . .« 72 
V. Severs Ovo eke RECOMMENDATIONS . « « s&s «© 9 « « ss 19 


BEPENDEX AS YWEANS AND STANDARD DEVIATIONS OF 
ecn Menge OARS BO 6-6 ww fs «teense «© “seas sv OD 


Meee DtxX  B: GEODETIC POSITION OF HORIZONTAL CONTROL 
STATIONS e ° es e ° ® e @ e e ° e e e es e ® s 91 


Mees DtxX Cy EMPIRICAL PROBABILITY DENSITY FUNCTION 
PLOTS e es t ° e es e e e e td e e e e se s e « 


\O 
Se) 


BIBLICGRAPHY e e e iL) ° es e e e e e ° e 6 e e e e e e e e 98 


MOPS ets ON LEST . 2. 5. 6 se ee lw tl ltl wl el ew!) CIO 





Oe 
Iti. 


nye 
Ve 
VI. 
722. 
wel. 
IX. 
Xe 
a. 
cit. 
Mert. 
XTV. 


Diol "Or, TABLES 


Sisciwiea esr ror Formulae —. < 


Data Acquisition Sequence of Evants ... 


Deaiee ngs rOr Standard Deviation (pooled 


SSaimaces) . 7% 


Summary of ANOVA Results at 95% Confidence 


Results of Interpolation Evaluation .. . 


Corrseted Aztrac Standard D2aviatio 
BExperinent Precision and Thasiolit 


Posz*ion Standards 


3000 Meter Arc. 
1900@Heter Arc . 
WiOOmMster Arc . 
700 Meter Are . 
500 Meter Arc . 
300 Meter Artec . 


Comparison 


58 
63 
65 
75 
72 
78 
85 
86 
a7, 
88 
89 
0 





LIST OF FIGURES 


He Mistaasson OL eRANnge= AZMAUTH Positioning 


Pistia OM OteAztaac Shore Station 


Pic NAG miEan SMe eng Ulit «os. « 


DrewuG@remtal Probab. lity Curve... . . 


Eccentric Range-Azimuth Gaometry. . 


Concentric Range-Azimuth Ssomatry . 


EeceiweGee  HEEOL COMPeENSAatlLon = < . ~. 


eaCn@oall DSC and dymec « 


Penamdeo-Azemuch Position . 


Vases On Of Gyywe Probability . « . 


Example of Angular Interpolation. . 


Sketch of the Survey Area 
Lines of Position Observed 


Uncertainty of an Observed 


Mees sO l.ae Lon Of Angle Only 


Yaeccpolation of Angles andi 
One-Dimensional Differences 


Two-Dinensicnal Difterenca 


Probability Density vs. Err 


Results Compared to dyms . 


Distance 
Batween M 
Retween 
2) ae e = a 


Results Compared to 90% Probability 


Probability Density Plot: 
Probability Density Plot: 
Probability Density Pict: 
Probability Density Plot: 
Probability Density Plot: 
Probability Density Plot: 


300 m arc 
500 m arc 
700 m arc 


1000 m arc 
1500 m arc 
3000 m are 


10 
Val 
19 
25 
Pe | 
28 
Ps, 
5)8 
22 
34 
38 
4 uy 
45 
43 
50 
50 
66 
67 
68 
74 
76 
ge 
25 
94 
35 
96 
a7. 








I. INTRODUCTION 


A. THE RANGE-AZIMOTH POSITIONING METHOD 


The fundamental purpose of a hydrographic survey is 
defined by Ingham (1974) as being t> "depict the relief of 
the seabed, including all features, natural and manmade, and 
tc indicate the nature of the seabed in a manner similar ta 
the tcpographic map of land areas." He goes on to describe 
two factors defining a single point on the seabed: 

eR eg poss iin. 0 * the point in the. 
orizonta lane in, for example, latitude 


and longitude, gtid co-ordinates or angles 
and distances from Known control points. 


(ii) The depth of the point below the sea 
surface, corrected for the vertical 
distance between the point of measurement 
and water level and for th2 height of the 
tide above the datum or reference level to 
which depths are to be related." 

Thus the hydrographer must answer the two primary ques- 
tions of "how deep" and "where" for each of the thousands cf 
scundings acquired cn every survey. Because every area tc 
ke surveyed has different gecphysical charact¢ristics and 
levels of use, the hydregrapher must possess a suite of 
tools and techniques to accomplish 2ach survey. A survey of 
a large metropolitan harbor requires different equipment and 
measurement precision than one for a deep ocean area. 

Only the first of Ingham's two factors cited abdcve is 
considered, and it is further narrowed in scope to *echni- 
gues used in the most precise surveys. Such a survey might 
be one of a winding, narrow river carrying deep draft 
vessels, or perhaps a very large scale survey of an inner 
harbor. Both areas require the highest positioning accuracy 


and 92 @inim@um of shore control stations. 





Any method of positioning employs the intersection of 
Mees Of position (LOP"'s) to construct a fix. Although 
advanced methods may use multiple LIP's, traditional hydrog- 
raphy uses the simple intersection of two lines to fix the 
vessel's position. The vessel is located somewhere along 
each of two lines of position, and the only point satisfying 
these conditions is the intersection of the lines. 

The error associated with one of the simplast posi- 
tioning methods, that of the two LOP range-azimuth fix, 
which is illustrated in Figure 1.1, is analyzed. 
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Plrguce 1.1 Illustration of Range- Azimuth Positioning. 


Also called the rho/theta method, range-azimuth positioning 
consists of the observation of a distance and an azimuth to 
a vessel from either one or two known locations (Umbach, 
1976]. An example of this method is the use of radar aboard 
Ship. A relative position for a radar contact is determined 
by observing 2 radar range and azimuth, or &@ tadar Tange ana 


weewal azimuth, to a contact. The two lines of position 


10 





always intersect at right angles because the observation is 
mad from a single point, and this concentric geometry 
provides the strongest fix possibls. Mariners also knew 
that the fix obtained via a visual azimuth is stronger than 
the one using a radar azimuth, because the visual bearing is 
more accurate. 

This example shows the advantages that make the range- 
azimuth method popular for hydrography. It provides the 
gecmetrically strongest possible fix, and only one location 
on shore need be occupied to control the survey. Such a 
positioning method is ideal in harbors or rivers where 
maximum accuracy is needed but where obstructions make cther 
types of fix geometry impractical. In 1982 the U.S. 
National Ocean Service (NOS) obtained twenty thousand linear 
nautical miles of launch hydrography, and sixty percent of 
this was controlled by the range-azinmuth method [{ Wallace, 
o's 3h. 

There are limitations associated with this method just 
as with any fix geometry. It is labor intensive and 
requires more radio communication (to establish fix timing) 
than mcst other methcds. In totally nonautomated situ- 
ations, distances to the survey vessel are recorded manually 
aboard the vessel, and azimuths ar2 recorded ashore by the 
theodolite observer at prescribed intervals. These fix data 
ace later put into ccmputer compatible digital form via a 
process called logging. Manual recording of these fix data 
are génerally toc slow to position every sounding. 
Therefore, individual sounding positions must be 
interpclated from the observed fixes. 

Systems have been designed that have an intermediate 
level of automation. The NOS Hydroplot System is an example 
of this type “Wallace, 1967]. When used in the range- 
azimuth mode, the vessel is usually steered along arcs of 


constant range from the theodolite station, and the 


ui 





Hydroplct System autcmatically records a distance measure- 
ment for each sounding. Azimuths are not telemetered to «he 
vessel but are relayed over voice radio and are manually 
entered into the computer system. Since the maximum data 
rate is akout two angles per minute, the interpolation of 
angles for scunding positions between fixes is necessary. 

Recently a digital theodolite, the Odom Aztrac, has been 
developed which can record and telemeter angles with great 
speed -- up td ten angles per second [Odom Offshore Surveys, 
Inc., 1982]. A computer system aboard the survey vessel can 
thus record and plot an observed position for each scunding. 
iaeamaaplad position fixing, combined with a computer's 
Seelity tO provide cross-track error indications to «he 
helmsman, enables the hydrographer to systematically cover a 
survey area with maximum efficiency by running straight and 
parallel sounding lines. 

The Aztrac system is still considéered a semiautomated 
system because an observer is required to manually track the 
vessel with the theodolite. Two fully automated range- 
aZimuth systems which feature fully automatic tracking have 
been devolcped. One is the Polarfix system developed by 
Krupp-Aclas Elektronik in Germany (Smith, 1983], and «he 
cther is the Artemis system develop2i by Christiaan 


Huygenslaboratorium in the Netherlands {Neweil, 1981]. 


B. HYDROGRAPHIC POSITION BRROR STANDARDS 


Historically, most national hydrographic organizations, 
qoweteeras the Enternational Hydrographic Organization 
(IEO), have used linear plotting error at the scale of tne 
survey to be the standard for sounding position accuracy. 
Prior to 1982, the standards recommended by IHO [IHO, 1968] 
were; 


Wimewindacated repeatabylity of a fix 
(accuracy of location referred to shore 


Ws 





ccntrol) in the operating area, whethe 
observed by visual or el@ctronic metho 
ccmbined with plotting error, shall se 
exceed 1.5 mo 50-05 in) at the scale o 
the survey." 


is 
asi. 
ldom 
fe 


The THO recently published new recommendations for e 
eeamea=aS [ITHO, 1982] which are: 


i, 


cor 


We... any probable error, measured relative 
tc shore control, shall seldom exceed 
twice the minimum plottable error at the 
scale of the survey (normally 1.0 mm on 
paper) .' 

Neither of the IHO standards make any reference as to 
what probability level they apply. Munson (1977) inter- 
preted the words "shail seldom exceed" in the above state- 
ments to mean "less than 10% of the tine", which seems 
reasonable. The 1982 IHO standard is somewhat confusing due 
to its use cf both the terms "seldom exceed" and "probable 
error", The latter term 1s associated with a 50% prob- 
ability by most statisticians including Greenwalt (1971). 
However, the author of these standards, Commodore A. H. 
Cooper, RAN (retd), has stated that he intended no statis- 
tical significance *tc the term "probable error" f{ Wallace, 
foe |. 

The NOS has not yet incorporat2d the latest IHO stan- 
dards, but such acticn is being considsréed in sone form 
{Weallace, 1983]. Current NOS standards have been developed 
to ensure that “accuracies attained for all hydrographic 
surveys conducted by NOS shall 2qual or exceed the specifi- 
caticns" c£ the 1968 IHO standards “Umbach, 1976]. Unlike 
the international standards, the NOS standards for all e¢leéc- 
tronic pesitioning systems use the concept of root mean 
Square error (dyms or rmse}), which has a somewhat variable 
probability of between 68.3 and 63.2 percent. fhe NOS stan- 
dards for fully visual and for hybrid (combination e¢lec- 
tronic and visual) pesitioning hav? no explicit reference to 


probability. 
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Specific 2perational standards for range-azimuth posi- 
tioning have been neglected by many hydrographic organiza- 
tions. However, NOS [Umbach, 1976] requires the following 
observational procedures be follow2i for all range-azimuth 
positions. 

"Objects sighted on should be at least 500 

m from the theodolite... the azimuth 

cneck Should not exceed one minute of 

arc... observed azimuths or directions to 

the sounding vessel for a position fix 

Shall be read to the nearest 1 min of arc 

Or better if necessary to produce a dosi- 

tional accuracy of 0.5 mm at the scale of 

che survey.‘ 
Since the range-azimuth method is classified as a hybrid 
positioning system, it is not referenced to any particular 
probability, but a reasonable assumption may be made that 
the dyms concept also applies in this case. 

The U.S. Naval Oceanographic Office (NAVOCEANO) also 
reguires that its surveys meet the standards of the NOS 
Hydrographic Manual. The Army Corps of Engineers presently 
have no formal positioning requirements that must be met by 
aii districts, aithough draft specifications are being 
Weeeem at chis «ime (Hart, 1983]. The range-azimuth tech- 
higue and its applicability to Corps of Engineers surveys is 
discussed in Hart (1977). Wo specific requirements for 
range-aZimuth positioning could be found fer either the 
Gaitemiman Hycrographic Service or the British Hydrographic 
Service. Palikaris (1983) also reports no published star- 


dards for these organizations. 


C. OBJECTIVES 


Alco sition Srror Standerds using an explicit prob- 
ability are based on the idea that an observation is a 
normally distributed random variable with zero mean and 


standard deviation go . These standards require a value for 


14 








Ea@emecraendard déviation of the component lines of positicn 
Pacem mare UL the £ix. Ths standard deviation is a value 
Suen hat there is a 68.27% probability of an observaticn 
falling within +190 cf the mean. It is unfortunate that 
often the 190 values of a hydrographic neasurement are 
Simply not known, or known only for ideal conditions and 
provided by manufacturers who have a vested interest in the 
measuring instrument. 

This paper and experiment, then, has as its orimary 
purpose the determination of a pointing error standard devi- 
atzon for two theodolites used to measure azimuths under 
hydr¢cgraphic corditicns. One instrument (the Wild T-2) is 
the standard used by NOS field units. The other (the Odon 
Azzrac) is a new digital telematering theodolite, which may 
prove useful in automating the presantly tedious manual 
methods. No 2Sstimation of standard deviation has ever been 
made on these devices under typical range-azimuth 
Genmdi tions. 

The T-2 1s the standard instrumant used by NOS surveying 
pazties for land surveying and, to a lesser extent, hydrog- 
raphy. It 1S a very precise instrument used in third-order 
horizontal control and provides an angular resolution of one 
Seecone of arc. However, it dces possess features that are 
less than ideal for range-azimuth work. The horizontal 
tangent screw of this instrument is awkward for range- 
azimuth surveying because it is not an infinit2? gearing 
device. The observer often encounters the end of the drive 
mechanism, stopping the instrument's movement while tracking 
the vessel. A soluticn to this problem for many observers 
is to track the vessel with the tangent screw unclamped, 
then clamp the screw cnly seconds before the fix cccurs. 

The inverted image feature of the T-2's in use by NOS 
Simplifies the optical system and reduces optical error. 
This is satisfactory for land survey work but creates some 
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confusion when tracking a fast moving vessel. This is 
because the obsetver sees an image of the survey boat upsids 
down and moving apparently in the opposite direction fron 
its actual movement. Another disadvantage is that the mech- 
anism for reading the horizontal circle of a T-2 is more 
complicated than desired for the rapid observations neces- 
Sazy in hydrographic survey work. The operator is required 
to step tracking the vessel, remove his eye from the tele- 
scope, and us? an auxiliary eyepiece to read she angle. 
There is a practical limit of about 30 seconds to the spesd 
Wich which successive angles can be observed and read. The 
pacticular unit tested was serial number 30504. 

The Odom Aztrac theodolite 1s a semi-automated, 
line-of-sight angle measuring system. The Aztrac system 
consists of a Wild T-16 *heodolite (Serial number 2534880 
was tested) which was modified for infinite tangent drive 
and te provide angular information in a digital format. The 
shore unit decodes the observed angle, determines the direc- 
tion of rotation of the instrument and displays the angie cn 
its front panel. The angle is then converted to binary 
coded decimal (BCD) format and used to frequency shift key 
(FSK) an FM transmitter to link the data with the survey 
vessel. The data is transmitted at the rate of 10 angles 
per second. On board the vessel the Aztrac receiver 
converts the received data to parallel form and displays it 
on the front panel for manual recording. For automated 
recording or processing by onboard computer a serial data 
oudut is provided. The Aztrac has an angular resolution of 
0.91 degree (36 are seconds) (Odom, 1983]. 

With the notable exception of angular resolution, this 
theodolite is more appropriate for range-azimuth hydrography 
than the Wild T-2. It has an erect image and infinite 
tangent screw which allow the vessel to be constantly 
tracked. Its digital output requires no action on the 
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Figure 1.2 Illustration of Aztrac Shore Station. 
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Operator's part to record observed angies. An additional 
positive feature of the digital output is the ability to 
tapidiy zero the instrument on an initial pointing. This is 
done by pressing a manual reset button on the instrument 
Gare to) panels) Figures 1.2 and 1.3 show photographs of the 
AZtrac equipment. Figure 1.2 illustrates a typical shore 
station. The observer is adjusting the Aztrac instrument, 
gieeecie Trans®atter unit is on the ground to the right of 
ene Agtrac trivod. The distance measuring equipment is 
mounted on a tripod behind the observer, and the Aztrac 
transmitting antenna is at the top of a pole on the extreme 
Tight of the photograph. Figure 1.3 shows the Aztrac trans- 
Beetee UNAt, With its digital angular display in hundredths 
of a degqrse. Both phctographs wer2 provided by Odon 
Omesmonre Surveys, Inc. 

A secend objective of this papar is to evaluats the 
interpolation methods used by the NOS for range-azimuth 
work. The a€vailability of the digital theodolite, with its 
direct measurament of all positions, enabled a comparison of 
observed and interpolated fixes to p23 made. This aade 
possible an estimate of whether interpolated positions meet 
required accuracy standards, and whather direct mnéeasureméent 
of all positions is needed. The estimate presented here was 
Mot Wade statistically rigorous so that the thesis could be 
kept to a manageable size. More theoretical statistical 
work is needed to fully reduce the int2rpolation data. 

The final objective of this investigation is t9 compare 
the position 2rrors of these two instruments with the 
Various position error standards discussed in section B of 
wrmomonapce:. The conclusions resulting from this objective 
Will assist the hydrogrephes to select equipment and oper- 
ating conditions that meet required position error 


standards. 
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Figure 1.3 The Aztrac Transaitting Unit. 
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II. ERROR INDICES AND RANGE-AZIMUTH GEOMETRY 

The method of range-azimuth positioning is usually 
selected for large scale surveys because of its simplicity 
and accuracy. This is the result of both angle and distance 
measurement devices being co-located, with the intersection 
of the two lines of position always being ninety degrees. 
In practice, however, the co-location of both instruments is 
of*~en not achieved. The result is an eccentric geometry for 
the position fix. This section will analyze the geometry of 
BPeemmeccentric and cencentric fixes. Position error indices 
in common use will be reviewed and analyzed for the special 
cases of range~azimuth methods, and th? error of an 


interpolated fix will be derived. 


A. DEFINITIONS 


Although a complete and general treatment of error 
theory will not be presented, some basic definitions are 
necessary to understand the data analysis presented here. 
The ideas in this section are included in many basic statis- 
tics textbooks, and were specifically drawn from Wonnacott 
(1977), Bowditch (1977), Heinzen (1977), Kaplan (1980), and 
Davis (1981). 

Error may be defined as "the difference between a 
specific value and the correct or standard value" (Bowditch, 
1977], or as "the difference between a given measurement and 
the "true" or "exact" value of the measured quantity" 


(Davis, 1981]. Mathematically it can be defined as: 


e= x,- T (2.1) 
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Mtere € 2S the esror, x;is an observation, and T is ths 
MCOLLTeEc~™’ Or "true" value. The word error ioplies that 
there is a known true value for a quantity, with which a 
measurement may be ccmpared to find th2 "error" associated 
with that measurement. Since the true value of a measured 
quantity is rarely knewn, the term "error" is not precisely 
correct. Davis (1981) states that it is more appropriate to 
speak of the theory of observations rather than che theory 
of errors, but it can be shown that the difference between 
the two is largely one of semantics. 

A Single neasurement of a particular quantity may be 
considered sufficient for many purposes, even if it is known 
that additional measurements will probably be slightly 
different than the first. If the quantity to be measured is 
of sufficient importance, then multiple measurements are 
made and the sample mean, X, is used. Each of these 
multiple measurements can be a considered numerical value 
for a random variable. A random variable is one that takes 
on a range of possible vaiues, each associated with a 
Bascmcular probablility. 

The sampl2 mean may be expressed mathematically by equa- 
Sone2.2 {acCnhnacott, 1977]: 


n 
xX = = ee (2.2) 
Asi 
where n is the sample size. If th2 sample size were 
increased without limit (n — oo), 2quation 2.2 would give 
the population mean «. The sample mean is always an esti- 
mate of the populaticn mean, which is never directly 
computed. This leads to the concept of the residual, v, 
which is the difference between the estimate X of the popu- 
lation mean and the observation x,. This is shown in 


Samat on. 2. 3. 
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ve X- x, (2.3) 


The residual is computationally the negative of the 
error. Nevertheless, equation 2.3 is more appropriate 
because it uses an estimate, X, of the unknowable pepulatic 
mean «a. The presence of X in equations 2.3 implies «hat 
multiple measurements have been mad2, and allows @ partic- 
ulér confidence +o be assigned to the estimate ae 
depending on the number of such measurements. Because <«=he 
word error is stili used in much of the hydrographic profes- 
sicn, it will be used interchangeably in this paper with the 
term residual. It is important, nowsver, to understand that 
the concept of the residual, whatever its name may be, is 
fundamental to any measurement oper2tion. 

Errors are classically divided into threes groups: blun- 
ders, systematic errer, and random 2rrors [Greenwalt, 1962}. 
Bowditch (1977) and Davis (1981) do not classify errors as 
wmemmeang blunders, but like the term error itself, the 
dietincticn is largely a semantic one. Ideally blunders and 
systematic errors are completely eliminated from the data. 
The most precise measurements reduce random error as much as 


possible, but it can never be completely eliminated. 


1. Blunders 


Blunders are mistakes that are “usually gross in 
Magnitude ccmpared to the other two types of errors" [{ Davis, 
1981], and ara most often caused by carelessness on the part 
of the observer, or by grossly malfunctioning observing 
equipment. They are usually detected and eliminated by 
procedural checks during the data acquistion process. The 
recognition of a blunder is not always easy, since a blunder 
‘may have any magnitude, and may be positive or negative" 
feemartcch, 1977). 
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2. Systematic Errors 


Systematic errors are defined by Davis (1981) as 
those that occur “according to a system which, if known, can 
always be expressed by mathematical formulation." This 
mathematical model results in correctors that are applied to 
all méasurements obtained, thus eliminating the systematic 
errors from the observations. The model may be as simple as 
a constant corrector subtracted from lengths obtained with a 
steel tape, or it may be as complicated as modelling the 
etfects of atmospheric refraction on electronic distance 
measuring equipment. 

If the systematic error is such that it cannot he 
modelled, it is then estimated by a process known as cali- 
bration. Kaplan (1980) defines calibration as the process 
of comparing the measuring instrument against a "knewn" 
Standard. The word "known" is usually operationally defined 
as a measurement Operation or instrument that is much more 
accurate than the one being calibrated. The difference 
between the observed and standard value is used as an esti- 
mecemomecnc cotal effect cf all systematic errors present. 
This precess is very clese to the classical concept of 
"errors" presented abceve, and is entirely proper for use in 
the correction of systematic errors { Davis, 1981}. Of 
course, one must be careful to apply the correcter cnly te 
those méeasurments made under the same conditions as the 
eae bration. 

A systematic error found in theodolit? cr sextant 
ebservations iS known as the personal error of the observer 
(Mueller, 1969], [Bowditch, 1977}. This type of error is 
rarely quantified for hydrogaphic applications, but never- 
theless it does exist. The observer must rely on the senses 
ef hearing and visior to make measurements, which vary 
between individuals as well as with time in one individual. 
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Some personal errors are constant and some are erratic 
(Davis, 1981]. These errors are minimized by training and 
standardizing observational procedures. The best way to 
eliminate personal error is by the use of completely 
automated observation equipment. 


"Random errors are chance arrors, unpredictable in 
Magnitude or sign", and are "governed by the laws of prob- 
ability" (Bowditch, 1977]. If one assumes that all blurders 
and systematic errors have been removed from the observa- 
tions, the remaining values can be regarded as sample values 
for a random variable. As noted earlier, a random variable 
can take on a range of values, each associated with a 
particular probability. A random arror has high probability 
of being close to the population mean, 4 , and a low preb- 
ability of being very much different than _« [Greenwalt, 
ioo2 ). 

A probability density function expresses the rela- 
tion between a value for a random variable and the prob- 
Sorwicy Of its occurrences. Hydrographic Survey measurements 
often use the normal or GausSian probability density func- 
tion. A concise 2xplanation of this function is given in 
Greenwalt (1962) and Kaplan (1980). The function itself is 
given aS equation 2.4, where p(v) is the probability of the 
eccurrence of a particular residual v, and g is the popula- 
tion variance which is approximated by the sample variance, 
S> given by squation 2.5. 
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p (v) (250) 


s= =) (Fx) (2.5) 
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Uemtanamenc probability of 2 residual falling 
between two residuals vy, and v., equation 2.4 must be inte- 
grated over the interval v, to v,. This corresponds to «he 
area under the gaussian curve between those two points, as 
is shewn in Figure 2.1. If p(v) ware integrated from -17 to 
+ig, the area under the curve would be 68.27% of the total 
area. This means that there is a 68.27% probability of a 
particular residual falling between plus or minus one stan- 
dard deviation of the méan, where the standard deviation,c, 
is defined as the square root of the variance given in 
equation, 2.5. 


==. ee See 6 See 2 a ee eee ee A SS SS Se EE ee 


=! 


p(v) 


‘ 


ee SE ey EE ge SE ee ey SE eee ie me ey lenin LN catia 


1 

oS owen aa c= 

<|—-—-—- 
op 

oO 

X< 

X 

ee 


: 
| 
| 
| 
| 
| 
| 
| 
| 
| 


Figure 2.1 The Normal Probability Curve. 


The combined effect of blunders, systematic errors 
and rand¢em errors can now be seen in overview. If it is 
assumed that blunders and systematic errors have been 
completely eliminated from a set of observations, there 
remain only random errors. If the sample size is large 


enough, then the sample mean and variance are good 
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approximaticns of the population mean and variance. That 
is, there is 2 68.27% probability that any future measure- 
ments made under the same conditions will not fall farther 
than plus or minus one standard deviation from the mean. 


B. TWO-DIMENSIONAL ERROR PIGURES 


The two-dimensional error of a position, as applied to 
the special case of range-azimuth fixes, must next be exan- 
ined. Two figures, the ellipse and the circle, are used to 
characterize two-dimensional error. The "error diamond" is 
also sometimes used but has no statistical significance 
(Thomson, 1977]. The error ellipse is discussed first since 
it is the most general index of error. Another is root néan 
Square distance (dyms), also known as root mean square error 
(Bowditch,1977], which is the radius of a circular figure 
commenly uséa in hydrceqraphy. It is the error index used 
for NOS positioning standards {Umbach, 1976]. A second 
Seweeteae Ligure, known as circular standard error, is also 
examined briefly because of the ease of converting it to 
circular figures which have different probabilities. 

Bee clarification of the issues involved in this 


section, the following assumptions are made. 


(1) Only random errors are considered. 
(21) Errors associated with each LOP are normally 
distributed. 


(iii) Errors are independent. 
(iv) Errors are limited to the two-dimensional case. 


(v) LOP's are straight lines at their point 


Gfe wunrersection. 
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These are the same assumptions made by Kaplan (1980), 
except that Kaplan allows systematic error to be considered 
in assumption (i). This appears to be an oversight since 
the remainder of his discussicn, from which this section 


draws heavily, considers only random error. 


Before proceeding further into a discussion of error 
figures, it is necessary to examine the two special cases of 
=e -azimuth pOSitioning, that of eccentric and concentric 


Geometry. Each is illustrated in Figures 2.2 and 2.3. 
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Figure 2.2 Eccentric Range-Azimuth Geometry. 


In actual practice the geometry used is often eccen- 
tric, but the concentric assumption is made. This is 
because it is usually difficult to co-locate both theodolite 
and ranging equipment. Hence, they are offset one or two 
meters from each other. It should be noted that this 
assumpticn will intreduce a systematic error in all 
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Figure 2.3 Concentric Range-Azimuth Geometry. 


positions, which must be 2liminated before an analysis of 
random errors can be made. The systematic error is often 
ignored (by use of the concentric assumption) if the total 
uncompensated error is within the tolerance of the standards 
being used. An algorithm for eliminating this error is 
shown in Figure 2.4 and given by equation 2.6. This algo- 
rithm assumes that the size of d in Figure 2.4 is small 


compared to fr. 


c= d cos 9 (2.6) 


where: 


the corrector to be appliad tor 


Q 
i) 


= observed range to the vessel 
distance between theodolite and ranging device 


©e@ a 
i) 


= the angle between the visual LOP and the 
line connecting the two stations 
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Figure 2.4 Eccentric Error Compensation. 
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Eeror Ellipse 


Detailed discussions of the development of the error 
ellipse can be found in many references, especially in 
Greenwalt (1962) and in Burt (1966). This paper will only 
present enough background to apply the error ellipse concept 
to wo LOP range-azimuth positioning. The error ellipse 
formed when multiple LOP observations are nade is not 
considered here. 

A range-azimuth position is formed by the intersec- 
tion of two LOP's, each having an associated standard devia- 
tion. By applying the two-dimensional normal distribution 
to the errors, elliptical contours of equal probability 
density are formed. The contours center on the intersection 
point of the lines of position. This is illustrated in 
Figure 2.5, and shown mathematically by 
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Ve-wGectatial (CLEOr) in ths disection oF on 
Weeoosmaawa lL) (eG595) in ths direc=-on of d, 
Jd. = langth of the semi-major axis 


x 
length of the semi-minor axis 
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M= Car be shown “hat 


2 A 
V. b 2 
—- + —7)=k (225) 


where 
K = - 2 tn | (va %) Cade 2 (2.9) 


Seeevalves Of pay ,v) from 0 to co , a family of equal 

probability density ellips2s are formai with axes K gd and 

Meese ctucs, '752), The probability density function in 

equation 2.7, when integrated over a particular area, 

Pacemes tne DS>bebility distribution function. This yields 
Ss 


q 


Bereprcbability that the r2 aiSev. anda View Wel "oeccut 
n 


d 
Simultaneously within that req nos OS Opabai it youd 2s - ta 


bee on sunc=2on of an ellipse is given by 


-Ky 
Pave.) = 1-— € (2. 10) 
Mieomsoru-20n of equation 2.8 for different values of K 
Melds d2ff€erant probabilities. For example, for 39.35% 


probability, the axes of the ellipse are 1.0000, and 1.000 
d, [Greenwalt, 1962]. In other words, a one-sigma error 
ellipse around a measured position indicates a 39.35% prob- 
Poe ieheveetnicew tne DOSItLON 1S actually within that ellioss. 
Tt is seen in Figure 2.5 that the standard devia- 
Pioweecmanasc, Of =he measured LOP*s are not the same as “he 


standard deviations d,and G of the error ellivss. A 
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Figure 2.5 Ertor Ellipse and d,,<. 
poesdoinae] transformation is =equiresd to obtain them. This 
SeancrOo>-mat"0n is found in Heinzen (1977) and will not he 


discussed here. Results of the transformation are presentad 
Beet auer2ons 2.11 and 2.12, from Bowditch (1977). 


g, = sag |tre 7 [(g, + a) Asie, (eis 
= | 2 2 ae 2 1 2 
Gns= sre |e e- (+ 0,)-4sin8e ;| (2. 12) 
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where: 
Ja = length of the semi-major axis 
Od, = length of the semi-minor axis 
% = standard deviaticn of the range LOP 
O02 = standard deviation of the angle LOP 


when converted to distance units 


es 


Pemerange-azwauth positioning, d, and dg, are not 


angle of intersection of LOP's 


equal. 0, is the error in distance measurement, and it is 
dependent on the @quipment used for ranging. A diagram of a 


range-azimuth position is shown in Figure 2.6. 
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Figure 2.6 A Range~Azimuth Position. 
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meemescOr in the visual LOP, d,, is a function of the 
angular srror of the theodolite and the distance to tho 
vessel, where r is the distance and Gis the angular stan- 
Gee d deviation "ee the theodolite, in units of degrees. This 
is given by equation 2.13, which is a modification fron 
Heinzen [1977]. Heinzen uses the tarm angular resolution in 
place of the more correct @&. 


r Jo 


Tis 57.296 


(2. 13) 

The error ellipse concept can now be applied +o the 
eccentric geometry by using equation 2.11, 2.12, and 2.13 
directly. In the concentric case, since the angle of inter- 
section (Be is always ninety degrees, equations 2.11 and 
2.12 can be simplified to equations 2.14 and 2.15. 
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Root mean square distance (d 1s presented here 


Pl 
because cf its common use in hydrography. It is not 
commonly known among hydrographers that unlik2 the error 
ellipse, dems has a variable probability depending on the 
eccentricity of its associated error ellipse, and ranges 
Promroc.3% to 63.2%, as shown in Figure 2.7 fron 

(ebatee 1977]. ECcentricity is defined as the ratio of the 
semi-mincr to the semi-major axes of the error ellipse. 
Root mean square distance is also called Mean Square 


Positional Error (MPSE) by Greenwalt (1962), who recommends 
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Figure 2.7 Variation of d,,. Probability. 


that this index not be used because of its variation in 
probability. 

Root mean square distance is defined in Bowditch 
(1977) as equation 2.16. An alternative form is given by 
Heinzen (1977) as equation 2.17. If the errors are assumed 
independent, the correlation coefficient Ga 1s zero and 
equation 2.17 is reduced to equation 2.18. An alternate 
method of arriving at equation 2.18 is to substitute equa- 


mmomec 2ekt and 2.12 directly into equation 2.16. 


ee 
Cie = 4 17, a oe (2ie: 16) 
| 2 2 
ad... = im 7 I,+ T+ 2a GG cos /3 (2.17) 
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dems * Sag Vt (2. 18) 


If concentric geometry is assumed, the angle 3 is equal to 
ninety degrees and equation 2.18 is reduced to equation 
ea 9 « 


dy ms -/6, eis (2. 19) 


Substituting equation 2.13 into equation 2.19, the final 
EOmpm fOr range~azimuth d,.,.. is obtained. 
2 is 


= g + 
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Circular standard error has experienced little use 
in hydrography but is valuable because it allows easy 
conversion between circles of different probability. It is 
derived in Greenwalt (1962) and given by equation 2.21. It 
should be noted that this equation is only an approximation, 
although a very good one. 


Gar=0,50000( 9, + D, ) (2.21) 


ircular standard error has a probability of 39.35% for a 
completely circular error ellipse. I+ is preferred over 
ieawemoecalse 1% can be converted t5 other circular error 
indices cf different probability by a constant conversion 
factor, as long as the ratio of %/dq is between 0.2 and 
imowmeethemequaticn fOr Circular standard error, and for 
other circular error figures, iS given in Table I, which is 
taken directly from Greenwalt (1962). This table gives all 


error indices in terms of either GG. or the error ellipse 
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| TABLE I | 
| Circular Error Formulae | 
| Percentage 
| Precision Index Probability Formula 
Circular 39. 390 = | 

| Standard Error | pe eee ( 0 + cy) | 
| WREN Gmjn/Omax > 0-2 | 
fe Cirouwlar 50% CFE = Geiltq iso | 
| Probable Error . | 
| GEE = 6.5887 (ee + oy) | 
| when oye Pre > 0.2 | 
| 

| CPE ~ (0.2141 Smin + 0-6621 max? | 
| when 0.1 < Geen ee Be ! 
| CPE ~ (0.0900 ore 0.6745 cn) | 

z 

| when 0.0 < 0, /o 1, < 0-1 | 
, Circular Map 90% CMAS = 2.1460 To 
Accuracy Standard | 
| CMAS = 1.0730 (a, + oy) 
| when ¢_._/o,4, > 0.2 | 
Circular Near- 99.78% 3.5000 og | 
| Certainty Error . 
: | 
| j 


semi-major and semi-minor axes Jd, and G,. Applying the 
assumptions of the previous section on dyms to this case, Q 
for the concentric range-azimuth cas? is given by equation 
Lie. 


a T 7 Ge 2.22 
f = 0.5000(%+ 55556 | ere) 
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C. THE ERROR OF AN INTERPOLATED FIX 


The Error in an interpolated position for both methods 
©£f interpolation used in range~azimuth positioning will now 
be derived. Fundamental to the derivation is an under- 
standing of the concept of error propagation, which is 
explained in detail in Greenwalt (1962) and Davis (1981). 
Error propagation is summarized here for the special case of 
range-azimuth positioning. It must be noted that any inter- 
polation discussed in this chapter is strictly due to errors 
in the observed positions between which the interpolation is 
made. Error due to the vessel not being at its interpolated 
position (due to steering or wind and séa conditions) will 


be considered in later chapters. 


1. Interpolation Algorithns 


The present NOS methods of interpolating range- 
azimuth fixes are of two types. One interpolates both the 
range and angle between two observed positions. The second 
is used when actual range information is acquizted on each 
sounding. In this case only the angle is interpolated, and 
is used with the observed range to compute a position. [In 
each case a linear interpolation is used [ Ehrhardt, 1979]. 
Algorithms for the interpolated value of the range, ©t,, and 
the angle 6;, are given in equations 2.23 and 2.24. 
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meted ita) (ood) (2. 23) 
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@, + (4/K+1) (9,-8, ) (2. 24) 


where; 
Subscript i denotes interpolated. 
Subscript 1 denotes observed position number 1. 
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Subscript 2 denotes observed position number 2. 
K denotes the number cf interpolations. 
j steps from one to K. 


The following example of the interpolation process 


is also sketched in figure 2.8. 
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Figure 2.8 Example of Angular Interpolation. 
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The values of the interpolated angles and rangés are 


computed by equations 2.24 and 2.23. 
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251 + (1/3) (244-251) = 251 - 2.3 = 248.7. (2. 25) 


Al 


251 - 4.6 = 246.4 (2. 26) 


251 + (2/3) (244-251) 


9; 2 
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1501 + (1/73) (1485-1501) 


1501-5.3 = 1495.7m (2.27) 
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ae 1501 + (2/3) (1485-1501) = 1501-10.7 = 1490.3wm(2. 28) 
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IEerOr Propagation 2s better tented the propagation 
of variances and covariances" [Davis, 1981]. The following 
paragraphs are a general derivation of error propagation, 
and will be applied below. In reading this section, the 
terms x andy are gener2l, yet ths reader should remember 
that they will be applied specifically to rangé-azinuth 
interpolaticn. Let y be a set of quantities each of which 
is a function of another set of random variables x. The 
random variables in cur application are the observed x 
Geom tn equation 2.1, and y is the interpolated ©, or fC. 
given in equations 2.24 or 2.23. 

The covariance matrix yey is given by the matrix 
equation 2.29, where Jy, is called the Jacobian matrix and 


is given in equation 2.30. 
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a 
oy. PSscne <ranspose Oe, and dixx iS another covariance 
Matrix given by : 


Yet Ors S21 GC, ! | (2537) 
l 1 
Orn, meats Ss ~ S,, 


The covariance matrix in eguation 2.31 is written 
for the general case cf correlated random variables. If y 
is a single quantity rather than a set of guantities, and if 
the variables are assumed uncorrelated, equation 2.29 


reduces to: 


0 oy 
2 Ox, 

on 
x aye 


els x. op 3X2 


dy Ox, ae x, C) we 
O, || | 


| 2X py 


cere y- fq Olt the multzplication in equation 2.32, the 


(2.32) 


expanded form is given in equation 2.33. 


2 


2 2 
2. p) = P) 2 d 2 
Jy = (%) 4 +(% C2 tr --+-- +> = Ce (Ze 33) 


If y is a linear function of random variables, the 
partial derivative terms in equation 2.33 become constants. 
Thus the matrix equation 2.29 becomes equation 2.34, where 


apomc, and d are constants. 


2 2 2 2 2 2 2 2 2 
Jy =agq + DO, +Coat---- + dq, izes) 
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Now let us apply equation 2.34 to the range-azinuth 
case. The assumptions made above are: 


(2) y is a Single quantity 


(ii) The random variables x are assumed 


uncorrelated. 
(arjeey ts a limear function of x 


ieee ness are applied to the two interpolation algcorithns, 


equations 2.24 and 2.23, it can be seen that: 
(i) @. and ©. are single quantities. 


(12) The random variables 9, , 98, and 


"=, . ©, are assumed uncorrelated. 


a2 


(lili) The two interpolation algorithms are linear. 


From the foregoing general iiscussion the interpo- 


lated variances of equations 2.23 and 2.24 can be given as: 


a 2 J > & a 
tA G,, + Cy (.- Tr) (2.35) 


Q 
l 





A yes 
Due, 
a 2 a 
Gan = Oe, = Gy (2. 37) 
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and: 


2 2 2 


Js, = %o2 = Te (2. 38) 


By SUbDStituting equations 2.37 and 2.38 into equations 2. 35 


and 2.36, we have: 


C= Oe (2. 39) 
and 
Ga = To (2. 40) 


This error propagation applies only to vessels 
moving in an arc. Fora vessel moving in a straight line 
the angular interpolation algorithm is not linear and the 
partial derivative terms analagous to those in equation 2.33 
are not constants. The error ellipse for an interpolated 
position can now be formed. The ellipse is ssen to be the 
same as fer an observed position (figure 2.5) since equa- 
tions 2.39 and 2.40 show that errors in interpolated 
distance and angle are thea same as for an observed distance 


and angl?. 
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III. EXPERIMENT DESIGN AND IMPLEMENTATION 

The experimental work performed in this investigation 
can be conceptually divided into two parts, although both 
parts were accomplished simultaneously. Part one involved 
pointing error and variance determination for the two theo- 
dolites, and oart two investigated the accuracy of 
interpolated positions. 


A. FIELD WORK 


The actual field work was reprasantative of a typical 
range~azimuth survey. Two full days (8 and 15 April, 1983) 
were required to obtain 443 position fixes and over 2500 
Aztrac angles. The experiment took place in southern 
Monterey Bay near the Monterey Harbor Coast Guard Pier, as 
shown in Figure 3.1. The vessel used was a chartered 
Bo-coct. Unt Elzte with a fiberglass hull and twin ergines, 
and its operator had about two months' hydrographic survey 
experience, including steering range arcs. 

For each position of the vessel, six lines of position 
were observed. Three Wild T-2 theodolites located at 
stations MUSSEL, SOFAR, ard USE MON were used t9 cbhtain the 
"reference" or best estimate positions of the vessel. Two 
additional test theodclites, a Wild T-2 and an Odom Aztrac, 
made observations frem stations T2 and AZTRAC. Finally, a 
Del Norte Trisponder (model ROY) provided a distance LOP to 
the vessel from staticn GEOCEIVER. The reference positions 
were obtained at one minute intervals, and five Aztrac 
angles were observed between each of these positicns for 
later use in evaluating interpolation methods. Figure 3.2 
shows the sources of the various LOP's to the vessel. It 
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Figure 3.1 Sketch of the Survey Area. 


represents an enlargement of the dashed circle area shown in 


Pegure 341. 


Geodetic control for the reference positions consisted 
entirely of monumented third-order stations. Station 
GEOCEIVER, although not a published station, was located to 
third order specifications by Mr. William Anderson of 
NAVOCEANO in 1982, and the other two stations were located 
as eccentrics of this station. Each was less than one meter 
from GEOCEIVER. GEOCEIVER is a monumented station, while 
Aztrac and T2 are marked by masonry nails driven into the 
concrete pier. Positions for Aztrac and T2 were computed on 
an HP-9815 conputer uSing the NGS geodetic ditect prtogran. 
The initial pointing by both instruments was to station USE 
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Figure 3.2 Lines of Position Observed to the Vessel. 


MON. Initial pointings at all stations were made to thizrd- 
Order staticns at least 500 meters distant. Geodetic posi- 
tions of all control stations is given in Appendix 8. 

Underway 2perations were very similar to a nonautomated 
Tangé-~azimuth survey, except that five theodolites, rather 
than one, were trained on the boat for each fix. The survey 
boat steered along the appropriate range arc, and £1x marks 
were given over voice radio to the observers on shore. All 
T-2 angles were obtained once each minute, and Aztrac angles 
were recorded every ten seconds. The sequence of events for 
each reference position is given in Table II. 

This process continued throughout the two-day field 
Operation. Breaks in data collection occurred at the ends 
of each range arc, and also when theodolite observers were 
rotated. The monotony of the events shown in Table II was 
sufficent to approximate the monotony (with its associated 
observer difficulties) of an actual survey. 
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| 1 
TABLE II | 
Slr: | 
Data Acquisition Sequence of Events | 
| ae ge Fe | 
i +o action { 
[ener ePrannnEl bar veesprareressuieererepseareveseemmer 
| 20 sec ‘7 depress Aztrac manual fix button | 
; { 
| 15 s¢c "15 second standby" over voice radio | 
| 10 sec | depress Aztrac manual fix but<+on 
| 5 sec ! "standby" over voice radio | 
: y 

Q sec "mark" over volce radio 
depress Aztras manual fix button, | 
T-2 Operators make observations | 
| | 
| | 
i 
a ee eran arense aap eareeny ane expen emeeee eneeneauns-cneameeen-enneenamnene expeuD aibeteanbene- anew) ammenseny cena eemensaarensans amass <D 


Observers at stations MUSSEL and USE MON were required 
to record angles as well as to observe. T-2 observers at 
staticns T2 and SOFAR were provided with separate recorders, 
because at times the high angular speeds of the vessel at 
these stations (due +o closeness of the boat) made observing 
mma Eecording difficult for one person. 

The Trisponder was calibrated using the standard NOS 
method over a gecdetic baseline (GEOCEIVER to USE MON) and 
ne systematic errors were observed. The geodetic baseline 
was determined by computing an inverse distance between 
Staticns USE MON and GEOCEIVER using the HP9815 computer and 
NGS geodetic software. The Trisponder was reported by Odom 
Offshore Surveys to have a standari deviation of its ranging 
error of 1.0 neters. The master unit and antenna were 
mounted at the highest part of the boat, about one meter 
above and aft of a radar antenna enclosure, which was about 
one meter in diameter and 0.4 meter high. The naster unit 
was a cube about 0.3 meter on a sid2 and was covered with 
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Green Signal cloth. This unit was the target to which 
angles were observed. The Trispondar Distance Measuring 
Unit (DMU) was mounted next to the steering station inside 
the boat. 

Data logging for both Trispond2r and Aztrac was accon- 
plished by an Odom Navtrace computer system. Both thea 
Trisponder DMU and Aztrac recéiver were interfaced to this 
unit. The Navtrace computer is programmed to automatically 
log fix data not on intervals of equal time, but on equal 
distance intervals from a reference line. For this reason, 
all fixes were legged by uSing the manual fix feature of the 
computer. This simply caused the computer to log distance 
and angle each time a button was manually pressed. Timing 
for fixes was provided by an NOS standard sounding clock. 
This is a mechanical clock with a buzzer set +o ring every 
ten seconds. 

The weather during both days of field work was good to 
fair for survey work. Visibility was good at all times, and 
winds were calm on each morning of operations. Afternoon 
northwest winds were a maximum of about 15 knots on bota 
days, which produced two to three f50t seas in the offshore 


part cf the operating area. 


B. THEODOLITE POINTING ERROR 


Meemoo-teing Srrot, x., -s defined as the difference 


between an observed and computed angle given by 
(3.1) 


and the mean pointing error, X, is therefore 


n 
pan | 
X = +) x; (37.2) 
A*1 
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where n is the number of observations, 6, is the observed 


position to the reference position of the vessel. The stan- 
dard deviation, G,, cf a set of pointing errors can be 
expressed mathematically as 


i Se ee 
par) (R-) 


Mats 2S the usual definition in most statistics texts and 


(3.3) 





agrees with the discussion in section II.A. of this paper. 
The ccmputed azimuth, 8, was dstermined by computing a 
geodetic forward azimuth from either test theodolite posi- 
tion to the reference position of the vessel. There is sone 
potential errsr in 6 due to an uncertainty in the reference 
position of the vessel. The reference position of ths 
vessel was determined by a least squares adjustment cf the 
three LOP's from the three Wild T-2 theodolites. A 
by-product of this adjustment is an error ellipse for each 
position. The angle subtending this error ellipse from each 
test theodolite is the error in 9@.. This is illustrated in 


Praguce=s. 3. 


jas 


: 


uncertainty 


A a a ca Ca, 
| Se ee a ED Ay Se, RS a 





Figure 3.3 Uncertainty of an Observed Error. 
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The least squares method was developed in 1794 by Gauss. 
It 1s based on the ccncept that as soon as redundant obser- 
vations are present there appear discrepancies among the 
observations. In order to eliminate these discrepancies a 
residual, v, nust be added to each observation. In the 
method cf least squares the residuals are determined such 
that the sum of their squares becomes a minimum, provided 
all observations have the same accuracy [Mueller, 1979]. It 
1s assumed that the three LOP's used in this adjustment do 
in fact have the same accuracy. Further information on «he 
method of least squares as applied to hydrography is given 
by Kaplan (1930). 

Reference positions were calculated by the least-squares 
FORTRAN program AZLSQ2, written by the author, which is a 
modification of program SILVA1 {Silva, 1979]. The progran 
generated plane and geodetic coordinates for each position 
fix, as well as the léngths cf the semi-major and semi-minor 
axes cf its error ellipse. The angl2 made by the semi-major 


axis and the x-axls was also computed. 


C. INTERPOLATION ALGORITHM EVALUATION 


The preseat NOS method of observing azimuths to the 
sounding vessel is totally nonautomated, as explained in 
Bec ven A. This precludes recording an azimuth for Each 
sounding because of the speed with which soundings are 
taken. Consequently interpolation algorithms are employed 
to plct soundings between observed positions. Iwo algo- 
rithms are used: one interpolates only the azimuth and uses 
an observed distance to the vessel, while the other intezpo- 
lates both azimuth and distance. Figures 3.4 and 3.5 


Syst Gpate the situation. 
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Figure 3.5 Interpolation of Angle and Distance. 


The method used to evaluate these interpolation methods 


fsesmmotle, Actual distances and angles to the vessel were 
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recorded every ten seconds, using the telemetering theo- 
dolite and associated computer logging system. A nenauto- 
mated system was then simulated by intarpolating between 
positions oktained every minute. Thus two data sets were 
produced, an interpolated set and a corresponding set of 
actually observed positions. These two sets of positions 
were then ccmpared and their differances examined. Both NOS 
interpolation algorithms were evaluated in this manner. An 


analysis cf the experimental results is given in Chapter IV. 


D. CHOICE OF EXPERIMENTAL CONDITIONS 


Any experiment must be performed under conditions 
Similar to those under which the results will be applied. 
This section will describe the specific error sources to be 
examined under range-azinuth conditions. By carefully 
choosing the experimental conditions, these error sources 
can be brought into the foreground for examinaticn, while 
all cther sources of error can be kept in the background. 

It is nécessary to intrcduce some standard statistical 
terminology to help understand the experiment design. The 
classical experimental metnod studias the effect of only one 
Variable. That is, it holds all effects but one to be 
constant, and varies that single on2 systematically. In the 
cas2= of this thesis, one factor (a theodolites) was varied 
systematically by introducing two lavels (Aztrac and T-2) of 
mene faetor. It is desired» to find the effect of that factor 
(the theodolite) on the error of a range-azimuth positicn. 
It is known, however, that the observer also has a consider- 
able effect on the error of a position. Thus another factor 

the observer) is introduced into the experiment. Four 
levels (four different individual observers) of this factor 


wer2e selected. 


a1 








The follewing four observers, who are attached «+o the 
Naval Postgraduate School, constitute the four levels of the 
ebserver factor. 


Peweeserald Eowttilis, NOAA. Instructor in 
Hydrography. Ten years expsrience in 
geodetic and hydrographic surveying. 


LT Maureen Kenny, NOAA. roi) eG plea 
hydrography. 2 years field experience 
atoard the NOAA Ship Davidson, primarily 
in Alaskan and West Coast waters. 


LT Mary C. Schomakar, NOAA. Student in 
hydrography. 4 years field experience 
aboard NOAA Ship Davidson, and on NGS 


horizontal control and leveling parties. 


Mr. James R. Cherry, Supervisory 
Geodesist, Naval Postgraduate School. 23 
years field experience in nydrography and 
geodesy, with th2 Naval Oceanographic 
Office. 


Although additional observers would have been desired, none 
wer2 available who had any experience with theodolites. 
Aside fron the factors to be investigated, lt was reécog- 
nized that there were background conditions shat affected 
the results of the experiment (Crow, 1955]. Some of these 
were taken into account explicitly in the design. The 
influence of all others was minimized by scheduling the 
experiment such that each combination of instrument and 
observer was evaluated in tcandom order, thereby randomizing 


the effects of these other conditions. 
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Blunders and some systematic errors were exp 
Seeounted for in the design. Blunders were iden 
= 


ally by their large magnitude and were simply delete 


(Qu 
th 
tf 
O 
=| 


the data set using the editor on th2 NPS computer. The 
Aztrac instrument essentially eliminates all blunders due to 
its automatic data logging feature. 

SyScematic errors integral to the theodolites, such as 
Somiemwerion eid eccentricity, were not accounted for, since 
they were very small compared to the size of the errers 
under investigation. One systematic error that could be 
Gjeee tatge is the initial pointing error, which is due 
largely to the observer. The NOS Hydrographic Manual 
jUmbech, 1976] requires the initial pointing =o be "accurate 
-ouwieenin £30 Seconds of arc". That is, the difference 
between beginning and ending pointings is not to exceed one 
Minute. For this experiment, the m2an of the beginning and 
ending pointings to the initial azimuth was algebraically 
subtracted from ail cbserved angles for that set. It could 
be argued that a more accuraté way of removing this systenm- 
aoac Error would be to prorate the difference in initial 
readings between the beginning and anding pointings. 
However, the former method is often used by NOS hydrographic 
Suevey units. 

One background ccndition, the angular speed of the 
vessel, was explicitly taken into account in the design by 
dividing the experiment into groups. Six separate subsets 
of the experiment were created by maasuring pointing error 
at six different angular speeds. These speeds were chosen 
+o be closely representative of speeds found under actual 
Smewevying conditions. To this end, integral vaiues of 
angular speed were not used, but rather integral distances 
from the vessel <9 the theodolite station. This corresponds 
to actual practice in most nonautomated or semiautomated 
Surveys. The vessel maintains constant engine speed, and is 
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navigated in a circle such that a constant distance value is 
displayed in the ranging equipment. For this experiment, 
distances of 300, 500, 700, 1000, 1500, 3000 meters were 
‘used. Using the approximate vessel speed of 2 meters per 
second (4 knots), the angular speeds in minutes per second 
corresponding to these distances were 22.9, 13.8, 9.8, 6.9, 
4.6, and 2.3, respectively. 

Morne: background condition that was explicitly 
accounted for in the design was the experience of the 
observers. For a completely general investigation, one 
would desire to evaluate the error using observers of widely 
varying experience. Only experienced observers were used in 
order to conserve resources of time and money, aithough ncene 
of them were experienced in using the new Aztrac instrument. 
All observers except LCDR Mills had acted as observer for 
range-aZimuth hydrography within the past year using the 
wma, T-2. 

There were other background conditions that affected the 
results of the experiment, including weather and lighting 
conditions, and observer fatigue. Other subtle factors may 
have also contributed to the error. In nonautomated 
systems, there might have been & time lag by the radio cper- 
ator in the vessel as the fix mark was relayed to the *heo- 
dolite operator ashore. The design of the experiment could 
not explicitly account for all these conditions, so several 
steps were taken to randomize the order of observations for 
each combination of instrument and observer. 

It was assumed that the theodolites used were randomly 
G@eawn from the entire population of T-2 and Aztrac instru- 
ments. This is a fairly good assumption for the T-2, since 
the particular unit tested is one of several maintained by 
the Naval Postgraduate School, and has seen séveral years of 
service. The assumption for the Aztrac theodolite 1s not as 
good, since it is one of only five in existence at this 
eigen Et was provided by the manufacturer for testing. 
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Theodclite observers ware required to be paired, since 
the experimental theodolites observed simultaneous angles +o 
the vessel. A random selection of obsarver pairs was accom- 
plished by consulting a random number table in Wonnacott 
V7 The Order Of observation for the six different 
range arcs waS important and required randomization. 

Failure to do so could allow the error associated with each 
to be influenced by time varying conditions. Examples of 
these are the changing effect of sun glare on the instru- 
ment, observer fatigue, the effect of repetition acting to 
decrease error, or increasing afternoon winds disturbing the 
theodclite. Randomization was again accomplished with the 
aid of a random number table. A random number was assigned 
+c each range arc, and the order of observation was estab- 
lished by selacting the arcs from the lowest to highest 


number. 


ays, 





A. DATA PROCESSING SYSTEM 


All data acquired for this experiment wer2 manually 
entered into the NPS computer facility. A data processing 
system was designed and implemented to log the raw data and 
perform the necessary computations to arrive at tne finished 
ferm given in this chapter. This system was divided into 
numerous subsystems to enter data, istermine the pointing 
error, comput? interpolated positions, compute means and 
varlances, test for randomness, and compute analysis of 
variance (ANOVA) statistics. ANOVA is a standard statis- 
tical technique used in testing for 2 difference between the 
errcetsS Of two Or BOTe Eactors. 

The data processing system was designed and programmed 
by the author. AS mentioned previously, the NGS gecdetic 
inverse and direct subroutines { Preifer, 1975], as well as 
the least-squares adjustment progran [Silva, 1979], were 
adapted for use here. Tha ANOVA computations were performed 
Dyeche Statistical Package for Social Sciences (SPSS) on the 
Naval Postgraiuate School mainframes computer {Hull, 1981}. 
The pointing error and interpolation subsystems were 


Gesigned es described in Chapter III. 


Be POINTING ERROR DETERMINATION 


The pointing errer standard deviation for sach theo- 
delite needs to be quantified, and a determination must be 
made as to wheather there is a statistically significant 
difference batween the two instruments. Pointing error for 
each combination of instrument and observer was determined 
by the methods discussed in Chapter III. Standard deviation 
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Gemeeie POINting Error was computed using equation 3.3, for 
each combination of instrument and observer, at each range 
arc, and are given in Appendix A. This entire procedure was 
performed first for the T-2 angles observed to the nearest 
second of arc, which is the maximum resolution for this 
instrument, and second, for T-2 angles rounded to the 
nearest minute of arc, which is specifically allowed by the 
NOS Hydrographic Manual [Umbach, 1976]. 

Since eight different instrument-observer combinations 
exist, the sample msans and standard deviations for each 
combination were slightly different because each sample is 
only an estimate of the mean, , and standard deviation, , 
for the entire population. The pooled standard deviation is 
the best estinate of for this case of multiple samples 
because it takes the differences among sample means into 
account. If an overall standard deviation is computed using 
equation 3.3, the population standard deviation will he 
overestimated because of the differences among sample means. 
The pooled standard deviation is mathematically expressed by 
equceson &.1 Crow, 1955], [ Box, 1978 ]. 





(ES <2 Grell 2 eee Ce ee 
(Aus) Shs Cos) Ss ie caliente De 


- (4.1) 
Nn, * A, +---- +n -K 


wheres 
k = total number of observer-instrument 
combinations 
h = number of observations for each obsérver- 
instrument combination 
S = Sample variance for each observer- 


instrument combination 


Thus a pooled standard deviation was computed for the eight 


samples available at each different range arc. These are 
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both the rounded and 


angular speeds for ths 


given in Table III for the Aztrac and 


unrounded T-2 data. This table gives 


vessel in units of minutes of arc per seconds of time. Note 


lets lpe 
that 
only 


the pointing error of the Aztrac is twice as large as 
See ric cmanametne: =ne rOoundea T=-2 poinving error is 


Slightly greater than the unrounded. 


| 
| 
| 
: 
| 
| 
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TABLE II! 
Pointing Error Standard Deviation (pooled estimates) 


| | 
| { 
( | 
| | 
{ Range Approx. ee, wee | 
| Are ee B2cEac (anrounded) (rounded) | 
ee 

(m) mae m (s2c) m (séc) m (sec) { 

| ( on Je { Ga do | Ga Js | 
ee || SS Pe |) ee rr 

| 300 | 22.8 | 3.33 (2290) | 1.26 (868) ! T.31.. °(9:0 2) | 
| 500 isso 35601470) 1.72 (712) | te 3 715) | 
700 | O28 | 22 Ome onl) 1.30 (382) 1.31 (386) | 
| 1000 6.9 2.97 (613) belie Sunt 2 5) | 1.99 (226) | 
| 1500 | bY (e | 4.25 (584) | 0.92 (126) | Oe 9357 ( 7128) | 
( 3000 | Ppa | 2-06 (142) | 1.41 ( 97) | 1.43 ( 99) | 
| | 
a a i es eee ween ee 


The technique of analysis of variance (ANOVA) 


is used 


determine whether a statistically significant difference 


exists between the two theodolites. 


The ANOVA technique 
Bivewsmene paz tlitloning of overall pointing error variance 


into portions caused by each factor (observer and cheo- 


Gol e) . 


by interaction, 


and by experimental error. 


Pieerfceronsexeses if the variance for a particular 


tion of instrument and observer is greater than the 


for any other such combination. 
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Paimanily due <0 he conditions described in Chanter III, 
Gemecerning the reference positions for the vessel. This 
error is a measure of the precision of the experiment, whére 
precision is defined as the closeness with which reneated 
measurements made under similar conditions are grouped 
together (Greanwalt, 1971]. Experimental error is assuned 
to have a population mean of zero. Further details of the 
ANOVA method can be found in many statistics textbooks, such 
as Wonnacott (1977), Box (1978), Crow (1955), and Walpol: 
(1978). This discussion was taken primarily from Crow 
1955). 

The following assumptions must be mad2 when using tha 


analysis of variance technique. 
(2) Observations are tandon. 


(li) Mears and variances ars additive, as given 


in the mathematical model below. 
(iii) Experimental errors are independent. 


(iv) Variances of the experimental errors 
are equal. 
(v) Distribution of the experimental errors 
1s normal. 
A mathematical model can now b2 given, using these 
assumptions, which specifies the total effects of the vari- 


ances acting on a particular observation, Xt: 


Riyt at az, +d, + 5:3 + @r-+ (4.2) 


4j 


wheres: 


“= overall mean for all observations 
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Lae cumceteoret he Ins -cument factor at level i 
by = effect of the observer factor at level j 


3; =NOtLece OF the Interaction of instrument and 


observer at level 1 and j, respectively. 


€é:.¢ = random effect caused by the variance of the 


experimental error. 


j-- 
i] 


1,2 
j= 1,2,3,4 
t = te 2, 3, @#@ 6 el 


n= the number of observations for a 
particular observer-instrument 


combination (usually 15) 


It should be emphasized thac <he a, b, §, and @ of equa- 
tion 4.2 ar2 not actual variances, but a realization of «he 
fect of those variances on a particular observation Xiit 

variances associated with observer, instrument, and 
interacticn are computed by thea ANOVA procedure, and these 
form the basis of the test for differances. This test is 
Gameg the F-test in honor of Sit Ronald A. Fisher 
[Wonnacot+, 1977], and is based on a ratio of variances. f[o 
test for a difference between instruments, we form a ratio 
with the variance amcng instruments in the numerator, and a 
denominator composed of an estimat2 of the variance of the 
experimental error. In térms of our mathematical model, a 
the numerator) is being compared with exit (the denomi- 
nator). More simply stated, it is a comparison of the 
precision of the instrument (the numarator), with the preci- 
SUOR Of the exoertment (the denominator). This ratio, F, is 
then compared to a ratio F, ,which is computed for for a 
particular confidence level from the Fedistribution function 
given as equation 4.3 [Crow, 1955}. 
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the ratio discussed above 


the number of degrees of freedon 


in the numerator of F 


the number of degrees of freedon 


in the denominator of F 


A precise hypothesis must now ba stated that can be 


tested by the F-test. Walpole (1974) states that "a statis- 


tical hypothesis is an assumption or statament, 


May not be true, concerning one or mor2 populations." 


Experiments are designed to test hypotheses, 
eeeon Or an hypothesis is to conclude that it is false, 
acceptance merely implies that 


otherwise" ff Walpole, 1978}. aA null hypothesis is an 


iweeetal hypctnesis, Or one we hope to reject" [Crow, 1955], 


and is usuaily stated in terms of an assumption of no 


difference between the effects to b2 investigated by the 


experiment. 


which are: 


(1) 


This experiment uses three null hypotheses, 


There is no difference betwsen observers, 
that is, the ratio F, for observers, is 


small compared to F,. 


There is no difference batween instruments, 
Cnate) “ge cnemrattour, Loe instruments, is 
small compared to F,. 


mere is No Pateraction, that is, the ratio 


Pil wre cacaclOn,s ls Stall Compared to F. 
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which may or 


and "the rejec- 
while 
we have no reason to believe 





These three null hypotheses can be accepted or rejected 
on the basis of the F-test and the ANOVA procedure. If any 
are to be rejected, the value of th2 ratio F must exceed 
that of the critical F, for that confidence level (95%). A 
Summary of ANOVA results 1S given in Table IV. For each 
range arc, the computed values of F and the correspcending 
critical value needed to reject the null hypothesis are 
given. The rightmost column of tha areas labelled as 
rounded or unrounded data indicate acceptance or rejection 
erecme null hypothesis. «It can be seen from this table that 
in no case could the null hypothesis be rejected for either 
instrument or observer. [In other words, it may be said with 
95% confidence that there is no reason to believe ther2 is a 
difference between the four observers or between the two 
instruments. It should be noted that data for the 1500 
meter are indicate a rejection of the null hypothesis for 
MeremeGr1On. ACCOfding to Crow (1955), a significant inter- 
action usually occurs because unrandomized background condi- 
sions are present. There is littl3 apparent reason why this 
interaction should occur, other than that there were some 
unidentified, time varying conditions affecting the measure- 
ments. The 1500 and the 1000 meter range are were only 
observed cn April 8, while all the other range arcs were 
investigated on both days of the experiment. 

A final assumpticn of the ANOVA technique is that the 
ObSerWations Within the eight conbinjations of instrument and 
cbserver are randomly drawn from thair populations. Several 
tests are available to determine if a particular sample is 
random. The test used here was the Run Test, as given in 
Crow (1955), and a confidence level of 95% was used. This 
test showed all eight samples were random at the 95% 
confidence level. - 
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TABLE IV 
Summary of ANOVA Results at 95% Confidence 


Rounded Data Unrounded Data 


| 

| 

| 

| 

| ee ote es | byp. 

| 50 so nn ee ana 
| ef | See | 8508 | core 
| observer | 20 69 ! 2522p" | Jc Oe accepe | 
| mmcecaction | 2.70 {| 41.36 | accept | 1.35 { accept | 
| Mico Secte sk ttti“(<i‘COSOSOSOSCO*~OCOOOOCOC~C—CC 
| re | | 0 | necoyt | 
observer 2.68 0.13 accept | 0.14 {| accept | 
| interaction | 2.68 | Ones. | AScepr | 0.25 | accept | 
eee ee! 
| instrument 3. &4 0:05 accept | 0.205 | accept | 
| observer {1 2.60 { 1.86 | accept | 1.86 | accept | 
| ieee reaction Ze OW 0.28 Sc S2De 0. 29 accept 
oe °° °° .°°.°.°.°.°.° ft 
! “instrument | 3.92 | 0.01 | accept | 0.01 | accept | 
| observer 2208) 1 UL 4o Mee ccept | 0.45 | accep 7 
| interaction | 2.68 1 1.74 | accept | 1.79 | accept | 
| 1500 m are | 
| “instrument | 3.84 | 0.78 | accept | 0.73 | accept | 
i observer { 2.68 { 1.44 | accept | 1.41 | accept | 
ineerac ion 2668 5.01 reject 5.02 | reject ! 
<i ||. 
| “instrument { 3.95 | 0.10 | accept { 0.12 | accept | 
| observer | Pape (3. | Uae, ascept 1.74 | accept | 
interaction | 2710 35 accept | le37 | accept | 
| | | | | 
| ea 
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C. INTERPOLATION EVALUATION 


The two methods of interpolation used by NOS were next 
evaluated. As explained in section III.C., three sets of 
positions wer2 computed. These were a set of actually 
observed positions using Aztrac, and two sets of corre- 
sponding positions computed by the two different interpola- 
tion algorithms discussed in section II.C. These two 
algorithms were evaluated by computing the distance between 
each observed pesiticn and each corresponding interpolated 
pose <20n. A FORTRAN program written by the author perforned 
+he cemputations. The NGS geodetic direct subroutine 
(Pfeifer, 1975] was used t9 comput? positions from cbserved 
distances and Aztrac directions to th2 vessel. The interpo- 
lated positions were computed using equations 2.23 and 2.24, 
which are the same algorithms used in the NOS interpolation 
Subroutine TCARC [{Ehrhardt, 1979]. Distances between the 
two corresponding positions were computed on a plane, rather 
than using a geodetic computation. There is negligible 
difference between plane and geodetic methods at the 


s 


ers. 


distances (about ten meters) under consideration : 
Results of these computations are given in Table V, which 
Shows for each range arc the average distance in neters 
separating the interpolated and the corresponding observed 
positions. An indication of the variability of these values 
is shown by the percentage of interpolated positions falling 
farther than 1.0 meter away from th2 actual position. 

Little can be inferred from the results in Table V 
because this table is really an intermediate step towards a 
rigorous evaluation of the raw data. This table should be 
viewed as only a general indication of the effectiveness of 
iMeemPelaticn. Since it is a direct result of the ability 
of the boat operator to steer the vessel in an arc, wind and 


séa conditions have a tremendous effect on interpolation 
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TABLE V 
Results of Interpolation Evaluation 


| 
| Range | Angle Only Angle and Distance : 
| | en > 1.0 m | a %> 1.0m | 
m X from mean {| X from mean | 
nn nn a nnn nnn nn nn wn nnn nnn ey 
| 3000 | 1.15 47% | 2.41 2% 
| 1500 fe 209 85% 4.98 91% | 
1000 ! EP loy>: 72% | 4.07 78% 
700 {| 2.04 75% } 3.01 80% 
| 500 | 1.86 19% | 4.15 3 3% | 
| 300 | 2.59 79% | 5.32 84% | 
| 
ee ee ee ee ee 


effectiveness. Althcugh this experiment was carried out in 
representative survey conditions, Table V should not be 
viewed as Deing applicable to all Situations. The table 
do2s indicate that, whenever possible, automatic recording 
of range data should be used. 

Full analysis of the interpolation algorithms should be 
the two-dimensional equivalent of tasting for the difference 
between means. This is because both the interpolated and 
cbserved positions are not "true" positions, bus have some 
error. A one-dimensional test of differences between means 
is well established, and is discuss2d in several references, 
including Wennacott (1977). The null hypothesis for such a 


test is 


d= A, - a (4.4) 
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where ma, and _m, are the means of tha two populations, and 4 
is some arbitrary distance selected by the experimenter. 

The two-dimensional problem has tha same null hypothesis but 
the mathematics of the test have not been established. This 
problem is illustrated in Figures 4.1 and 4.2. 


Q. 


SS ee ee ee eee 


| 


Figure 4.1 One-Dimensional Difference Between Means. 


A proper analysis of the data would inguires for 
interpolated-9bserved pair of positions, whether the 
distance between the two positions was greater than d for a 
Baeeeculat confidence. More work than could be incorporated 


into this thesis is required to fully evaluate the data. 


D. ANALYSIS OF FACTORS AFFECTING THE RESULTS 


The results of this experiment wera presented in Tables 
III, IV, and Ve. An attempt will now be made to analyze the 
experiment for errors in logic and technique, in order to 
better understand these findings. 
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Figure 4.2 Two-Dimensional Difference Between Means. 


The most lmportant results from this thesis are the 
estimates of g@ for theodolite pointing error. It is obvious 
from Table III that standard deviations of Deri ene «round ed 
and unrounded T-2 error values are about one-half that of 
the Aztrac, for all angular speeds considered. The ANOVA 
technique, however, shows no statistical differance between 
the instruments at the 95% confidence level. An analysis of 
the data used to obtain the results yields a potential 
Saelemation f25 this apparent contradiction. 

The original data (pointing errors in seconds of arc) 
were made the subject of empirical probability density plots 
using the subroutine HISTG [Robinson, 1974] on the NPS 
computer. These plots show probability density versus 
error, as well as mean and standard deviation. An example 
of these prlets, for the 500 meter range arc, is given in 
Figumes tas. The remagning plots ere found in Appendix C. 
Paets are shown tor beth Aztrac and,T-2 (unrounded), for 


each range arc. A striking featur2a of the Aztrac curves is 
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Figure 4.3 Probability Density vs. Error. 


their bimodal shape, as compared to &@ single peak for the 
T-2 curves. [+t can be easily seen from the curves how the 
spread of this bimodal distribution would increase the 
computed standard deviation for tha Aztrac data. 

The method of data acquistion for this thesis was 
semiautemated in that all T-2 angles were manually record:d, 
While the Aztrac angles were recorded by presSing a button 
aboard the vessel. It is probable that a time lag existed 
between all the T-2 observations and the Aztrac observa- 
tions, despite the best efforts of the observers, because 
the observation procedures was not totally automated. If 
this were true, there would be little difference between the 
observed angle © for the T-2 and the computed angle 9, 
because @, is associated with a reference position also 
derived from I-2 observations. The observed angle 8 fer 
the Aztrac would, however, be consistently different from 8% 
because of this time lag and because the vessel was moving 
to the left or right with respect «> the Aztrac observer. 
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Since approximately equal numbers of observations were made 
Meceethe NOd= moving to the left or right along «he range 
arc, the distribution of Aztrac pointing errors would take 
On a bimodal shape. 

Iz is understood that the data analyzed by these curves 
come from not one but four different samples, and thaz the 
curves should not be expected to ba perfectly peaked. The 
T-2 curves, however, also come from four samples and do not 
have multiple peaks. This analysis is further supported by 
finding the distance kLeéetween one peak of the Aztrac curve 
and the single peak of the T-2 curve in figure 4.3. fTh32 
distance in arc seconds, when converted +o meters, is 
roughly the distance the vessel traveled in one second. One 
second of time 1S certainly a reasonable figure for the tine 
lag discussed above. A manual check of the raw data 
recorded in the fleld also Suggests such a time lag. The 
original data were sorted into two sets of "left" and 
NMeight" cbservations, which were analyzed for mean and stan- 
dard deviation. Results of the analysis are shown in Table 
VIe This table gives the mean and standard deviation for 
the "left" and "right" data sets, and shows chat the mean of 
both sets was about two meters to the left or right of the 
reference pesition of the vessel. This two aster difference 
corresponds ciosely to a nominal vessel speed of two meters 
per second (four knots) for the boat used, and a time lag of 
one seccnd. Means for tne 1500 and 3000 meter rtanges are 
somewhat unequal because sea conditions at these offshore 
ranges caused «he boat to travel slower in one direction. 
The rightmost column in Table VI gives the pooled standard 
deviation of each "left" and "right" data set, which is the 
best estimate of the population standard deviation , 9, , for 


tHe -Azcrac. 
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TABLE VI | 

Corrected Aztrac Standard Deviation | 
| 

| 


| 
| 
( 
| 
| 
' | | Mean, X | Ss jAztrac Pooled 
| Arc = fon | Standard 
| , lr i ont m (sec) | (sec) Deviation | 
econ le 2.95 20:35 )a lel 7 806) | ( 

i 300 | | i ( y 1.07 (741) ( 
| Rk 1) =—2. 67 (-W836) | 0.99 (682) | 
| | o& {| 2.40 ( 993) {| 1.50 (660) | ‘ 
| 500 | 1.44 (598) | 
| i R jy 72.69 (71112) | 1231 (539) | 
( i Po 1.61 ( 476) { 0.90 (266) | ( 
| 700 | 120.5:° (312) | 
| i L {| 2.49 ( 514) 4 1.55 (321) | | 
i 1000 ( } 1.43 (296) | 
| R | -2.92 (-604) { 1.30 (268) | 
| iL | 1.93 (266) [io55 (2285 | | { 

1500 | } 1.45 (199) 
| | oR = -—-3.22 (-4 4) ( 1.096 (147) | | 

L 0.43 (367) ( 1.42 (98) | 

3000 | | | { 1.30 (90) 

| R ~3.97 (-273) | “Wee22 CSR | 


eel 


The systematic error caused by a time lag as the boat 
moved left or right in the observer's field of view was the 
result of fauity design of the experiment. The proper way 
+o correct this problem would be to duplicate the experiment 
using better synchronization of all observations. An alter- 
Rative would be to9 model the systematic error and apply 
corrections to. the existing data. Such a model should 
include an estimate cf the boat spead and its left or right 
direction with respect to the observer. 

Different ANOVA results might b2 obtained using the data 
corrected for systematic error, but this would not explain 
the ANOVA results in Table IV. Tha conditions affecting 


data acquistion must again be-considered, as well as an 
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understanding of the ANOVA process, when offering an 
explanation. 

The denominator of the F-ratio discussed in section 3B is 
essentially the experimental error of the measurement 
Bsecess. This error is primarily a function of the error 
ellipse for the "reference" position shown in Figure 3.3. 
Table VII gives an estimation of this experimental error, by 
Comparing the size of the error in 3, and 8.. Uncertainty 
in ®, is given as the méan major axis of all error ellipses 
for a given range arc. Uncertainty in 8, is the pooled 
Stan@ard deviation , G52 Ox each test eHeSOdGoL@ce. The 
values for 6, in the aAztrac column are from Table VI, and the 
H—2)Velues are from Table III. Am 2kXamination of Figure 3.3 
shows that use of the major axis of the ellipse is a worst 
case estimate of the uncertainty, since the ellipse could 
have any orientation in the x-y plane. Thus it can be seen 
from Table VII that the uncertainty in 8 1S smaller than 
that of the observed azimuth 6,. This comparison is an 
indicator of the precision of the experiment. 

If this error in the computed angle, &, could be 
reduced by decreasing the size of the error ellipse, che 
denominator of the FP-ratio would be smaller and the ratio 
itself weuld be larger. Thus a nore precise experiment 
could produce a rejection of the null hypothesis, altkougna 
thas is net indicated in iight of the values of dp, for the 
Aztrac and T-2 shown in Table VII. Reducing experimental 
error any further than this study would be difficult under 
typical hydrographic conditions, because the three LOP 
theodolite intersecticn position, adjusted by the least- 
squares method, is one of the most precise positioning means 
available tcday. The size of the reference position error 
ellipse could possibly have been reduced if better intersec- 
tion angles were available, but this was not possible with 
the particular geographic shape of Monterey Bay. 
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TABLE VII 
Experiment Precision and Iheodolite Error 











| 
| | 
| | 
Range, | Mean major RZemac pocled | T-2epeoled 
| meters axis std.deviation | std. dév. 
| meters (ec) meters (sec) | meters (séc) | 
ee 
| 300 ~974 (668) 1.07 (741) 1.26 (868) | 
| 500 -838 (344) 1.44 (598) ev 2 (€7172) { 
i 700 | 754) (222) Nie S55 C3 12) es'0ne( 382) 
| 1000 | ~824 (170) | 1.43 (296) 1.08 (225) | 

| 
WoO | . 738 (102) j We 5) (199) | 0.92 (126) | 
3000 | 1.2 SaaS 2) | 1.390 (90) | 1.47 (97) 
| | 
Teen gi et 


In summary, 
data’ do, 
ence between Aztrac and T-2 at 95% 


*he analysis of variance performed on these 
Piece, show that there le no Statistical dif fer- 
confidence. This proce- 
durée is a strictly numerical one and must be viewed in light 
of the original reasearch conditicns. When these conditions 
acre carefully analyzed to remove as much systematic error as 
possibl2, the data strongly suggest chat the ANOVA results 
are indeed correct. 


Ee APPLICATION TO POSITION ERROR STANDARDS 


From Chapter I it can be seen that there is some confu- 
Sion in the hydrographic community as to the application of 
With this 
there appear to be four possibilities for consideration as 


Pecpepttity tS positioning standards. abel fblse\ tal 


tandards with which to compare tha results of this thesis, 


as given below. 


(1) The 1968 IHO standard of 1.5 mm at the 


the 





Survey scale, with the 90% probability 
Suggested by Munson (1977). 


(11) The 1982 IHO standard of 1.0 mm at the 


scale of the survey, with 90% probability. 


(111i) The current range-azimuth standazd of NOS 
(0.5 mm at the survey scale) assuming @ 


DrObao Mt yeassOCiatecdmeitnh dak « 


(iv) The dyn; Standard of microwave range- 
range positioning found in Umbach (1976). 
This requires that dyms values at the 
Survey scale not exceed 0.5 mm for 1:20,000 
scale surveys and smaller, 1.0 om for 
1:10,000 scale surveys, and 1.5 mm for 


Surveys of 1:5,000 scale and larger. 


These four possibilities may now be compared to the 
Boomer 20N GCrErors of the AzZtrac and I[-2 by uSing the pointing 
error results given in Table VII. Since the precision of a 
range~azimuth positicn is also a function of the standard 
deviation of ranging error, 9 , assumed values cf 1.0 and 
3.0 meters are used in this aralysis. Some manufacturers of 
microwave ranging equipment used in range-azimu*h posi- 
tioning report a 1.0 meter value for oo, but 3.0 meters is 
most often used by NOS personnel (Wallace, 1983] and has 
some supporting experimental evidence by unbiased 
experimenters (Munson, 1977 }. 

Using these assumed values for gd, and the observed 
values for 9d, given in Table VII, the error ellipse axes may 
be computed using equations 2.11 and 2.12. However, the 
error ellipse must be converted to a circular error figure 
in order to be compared with the standards above. The 
assumed standards (i) and (i1) must use a-90% probability 
Circle, which can be computed using Table I, and standards 
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Figure 4.4 Results Compared to d,,p,;. 


(iii) and (iv) must use the d,,, formula given by equation 
2.16. Results of these computations are shown graphically 
fr oe 425 for 90% probability, and tm Figure 4.4 for 
dym; - Each figure has angular speed of the theodolite along 
the abscissa, and an ordinate consisting of a distance 


Scave, im @ec2rs, indicating the radius of the error circle. 
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[it] SeEror circle for Figure 4.5 is associated with the 90% 
Paopemality level, while the dyms circle for Figur2 4.4 is 
or somewhat variable probability, as discussed in Chapter 
II. The points plotted in Figure 4.5 have eccentricity 
Walues (92/, ) ranging from 0.3 to 9.6, which indicate a 
probability rang2 of from about 65% to 67.5%. 

These figures also clearly show an improvement in the 
estimate of Jg for the Aztrac as a result of the time lag 
correction discussed earlier. The position error values for 
the uncorrected Aztrac error is shown by triangles in both 
figures, while position error computed using the corrected 
Jo Values are shown by solid dots. For both Aztzrac casss, 
two solid linéar regression lines are drawn. Both figures 
show that the uncorrected Aztrac values have a much greater 
variability than the corrected values, and that the 
corrected values are almost the sama as the T-2 position 
error values indicated by open dots and dashed linear 
regressicn lines. If the corrected Aztrac values are taken 
to be the best estimate of position error for this instru- 
ment, then a relatively constant error is indicated for the 
entire range 9£ angular speeds considered here. This 1s 
about 3.3 meters dvyms and 4.6 meters (90% probability) for 
both instruments using a G value of 3.0 meters, and 1.6 
meters d,,,; and 2.5 meters (90%) for a G, value of 1.0 
meters. FPlots are not shown for computations using a range 
Senormor 1.0 Ret Sr. 

The four possible assumptions for position error stan- 
dards are compared to Aztrac and T-2 position errors in 
Table VIII. The roman numerals heajiing the columns of this 
table refer to9 the position standards associated with the 
Same numerals at the beginning of this section. The reader 
should use the table by selecting one of these columns and 
imepecting 1t trom the top to the bottom of the table. The 
first row of Table VIII indicates the probability associated 
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Figure 4.5 Results Compared to 90% Probability. 


With each positioning standard. The second row lists the 
maximum position error allowed by that standard, at the 
scale of the survey. Rows three ani four show the errors 
allowed in row two, when converted to actual distances for 
two representative survey scales of 1:5,000 and 1:1710,000. 


Rows five and six shew the radius of the associated 
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Peemoability circle fcr both Aztrac and T-2, for ranging 
error values 5f 3.0 and 1.0 meters, respectively. The dual 
probability percentages in columns (iii) and (iv) indicate 
the variable probability of d,.,. Dual values in cclumn 
({i1) for maximum error at the surv2y scale result from the 
NOS standard for range-range positioning. The remainder of 
the table presents conclusions as to whether the T-2 and 
Aztrac meet the varicus standards. For example, in column 
(iv) the observed 3.3 meter d,,, value in row five is less 
than the maxinum allowable error of 5.0 meters shown in row 
Bel@ene Tneretore, the T-2 and Aztras do meet the NOS range- 
azimuth standards of 0.5 mm at the scale of the survey for 


1:10,000 scale surveys. 


ai, 





TABLE VIII 
Position Standards Comparison 


| Fi | 35) u( 253.) { (iv) 

HO ( HO ( OS NOS 
emer ue CO Cem (te) Gey 
Assumed _ | | 90% {| 90% | 68%- | 68%- 
Prebability | | i 763% | £=-63% 
Allowable { | | \ 
am EXDror | 1.5 mom 1.0 “mm ( 1.5, ( 0.5 ann 
at Scale / 1.0 mon | 
Allowable | { { | 
Max Error room 8) 550 @ | 7.5 ol 2.5 8 
ae 125,000 | ( ( 

Allowable | ( ( 
Max Error 15.0 o | 10.0 m { 10.0 o 5 ol 
at 1: 10,000 | 
Aztrac & T-2{| | | 
EacOor i 4.6m 4.6 m { 3.3 o Sooo mt 
Meee! et el 
azecac & T-2{ | 
mer or 2-5 m | Zen m4 1.6 m {| 1.6m 
(oO = 1.0m) | { | 


@eoweeeatP eae eB se BSP BPP 2 SF Pe FSF SP SP BPs ses SF Fe -—s FBP seessei ates Fae SFB SS Bee ee ee ee ee ee Se Se 


mene eet) Blam (mE) 
Vso 000 | yes { yes { yes | no 
scale | | | 
1: 10,000 | yes yes { yes , yes 
scale i 


“4:5,000 | yes { yes | yes | yes 
scale 
1:10,000 yes yes 
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V. CONCLUSIONS AND RECOMMENDATIONS 

Of the original cbjectives for this thesis discussed in 
Chapter I, the first and most basic was the determination of 
Pelmeing error standard deviation for the Aztrac and T-2 
theodolites. No investigation of this type had ever been 
done for conditions typical of a range-azimuth survey. An 
experiment was carefully designed t9 determine this pointing 
error and to determine if there was a statistically signifi- 
cant difference between the instruments. The initial esti- 
mates of pointing error were given in Table III, which shows 
the Aztrac to have an error, when converted to distance, of 
about 3.0 meters, while the estimate was about 1.3 meters 
roan T—2 . 

An uncompensated systematic error in the data, due to 
the time lag discussed in Chapter III, was discovered when 
empirical probability density function plots were made of 
the entire data set for each instrument. This led toa 
revised estimate of the pointing error for the Aztrac, 
because the bimodal distribution causad by this time lag 
adversely affected the original estimate for g@. This 
revised estimate is about 1.3 meters, as shown in Table VII, 
and is almost the same as the value of G for the T-2 
instrument. When these estimates for dg are viewed in light 
of the precision of the experiment, as shown in the table, 
it is seen that the actual values of 0, could be smaller 
than indicated, because smaller values would be masked by 
the relative imprecision of the experiment. It is to be 
concluded, however, that the actual values of the pointing 
error of each instrument are no larger than those given 
here. 


Ee 





The question of a statistically significant difference 
between the instruments was then considered using the ANOVA 
technique, which can be said to compare a variance component 
due to the instruments with a variance component due to tke 
precision of the experiment. This precision was not very 
Much greater than the variance of the instruments, but was 
based on the most precise positioning method generally 
available for hydrography -- an intersection position using 
three theodolites. The ANOVA proceijiure indicat2d neo signif- 
icant difference between instrumenzs, put if the precision 
of the experiment had been increased, a significant differ- 
ence could possibly have been detected. In light of the 
subsequent discovery and elimination of the systematic error 
due to a time lag, this conclusion of no difference between 
theodolites appears to be well justified. 

An evaluation of interpolation aethods was the second 
objective of this thesis, and although the analysis was not 
as rigorous as it could be, it can be concluded that there 


is @ measurable distance between an interpolat2di position 


4) 


and a correspsnding observed position. This has never been 
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done for the case of range-azimuth positioning because ¢ 
rapid position fixing available with Aztrac has not been 
available. It has been shown, through an é@rror propagation 
analysis of the interpolation algorithms, that the interno 
lated error is not inherently due to the algorithms then- 
selves. The error is therefore due to the inability cr the 
vessel tc follow the range arc, whith is caused primarily 
because of environmental conditions and the vessel opera- 
tor's track keeping capability. The distance between «he 
interpolated and corresponding observed positions may be as 
much as two to four meters, as indicat2d in Table V, and is 
roughly twice as great for a position *hat is computed using 
both distance and angular interpolation, as for a position 
using angle interpolation alone. It is therefore 
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recommended that, whenever possible, automatic recording of 
range data should be used. 

The third and most important objective of the thesis is 
a comparison 9f the tctal position error using these instru- 
ments with the required error standards of the major hydro- 
Graphic survsy organizations. The lower half of Table VIII 
gives these results, which are that all the standards 
considered are indeed met, except the NOS range-azimuth 
Standard at 1:5,000 scale using a rang2 error of 3.0 meters. 
Bivs conclusion requires a very important qualification 
regarding the T-2 instrument and th2 arrors encountered 
While pursuing the first two objectivas. These are errors 
due to interpolation, and to the tine lag discussed in 
Grapter 19% 

The approximately one second time lag discovered with 
the Aztzrac data set is not actually associated with the 
Aztrac at all but is associated with the T-2. I+ appeared 
to be a systematic error of the Aztrac in this experiment 
only because the referenc2 postions were obtained using T-2 
instruments Similar to the test T-2. It must be concluded 
S=on the data acquired in this project that there exists, 
for any angle measured with a T-2, a time lag of about one 
second between angle observations and any measurement made 
aboard the vessel, including both automatic and manually 
recorded depth and range data. There is then an associated 
position error for these measurements, the magnitude of 
which depends upon vessel speed, which was about two meters 
for the four knot speed used in this experiment. The 
conclusions and position accuracies for observed T-2 posi- 
tions in Table VIII do not take this additional error into 
account. When the error contributions from both the time 
lag and interpolation are considered, it can be concluded 
that positions interpolated between observed T-2 positions 


have an additional error of about two to four meters. Thus 
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eae —<e6tal accual position error for T-2 positions and posi- 
tions interpolated between T-2 positions mignt fail to meet 
more of the standards than are indicated in Table VIII. 
Having reached conclusions related to the objectives of 
the thesis, another set of conclusions and recommendations 
can be made for the Odom Aztrac theodolite, regarding its 
ease of use and suitability for range-azinuth hydrography. 
The Aztrac instrument was expressly designed for range- 
azimuth or azimuth-azimuth positioning, and has feature 


Ss 
which are advantageous to the operator of the instrument 


rw 
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the field. Such advantages are discussed in Chapter I, and 
include ease of tracking, because of an upright telescope 
image and an infinitely geared tangan* screw. 

The Aztrac theodclite possesses advantages much more 
important than ease cf use in the field, and «hese addi- 
tional advantages are derived primarily from its ability to 
be interfaced with a computer aboard the survey vessel. 
Besides eliminating systematic error due to a time lag, a 
computer pased survey system offers the additional advantage 
of being able to measure and récord a position every few 
seconds. This allows three important advantages over a 
system that can only measure positions once per minute. 
First, since each position is individually measured, no 
interpolation is required and thus better accuracy is 
obtained than with a nonautomated system. Second, no manual 
data logging is required, which reduces blunders and greatly 
increases the speed with which a survey may be processed. 
Third, an automated system allows the surveyor to run 
straight sounding lines rather than curved arcs, because an 
automated system can provide an almost real-time cross track 
error indication to the helmsman. Running straight sounding 
lines increases survey efficiency by orienting the lines 
more normal to the depth contours, and by requiring fewer 
total linear miles of hydrography for each survey. A vessel 


82 





may of course steer straight lines with a nonautomated 
Syeorcm) DUtEAt 1S extremely difficult to maintain the stzric+ 
line spacing requirements for hydrographic surveys without 
an indication of cross-track error, so curved range arcs are 
usually fcllowed. 

Disadvantages encountered during this experiment incluée 
additional operator fatigue caused by the requirement to 
constantiy track the vessel. This is necessary if the 
substantially increased data rate available with this 
instrument is to be utilized. The extra effort to track the 
vessel is mcre than offset, however, by the elimination of 
manual data logging. Care is required by the operator when 
rotating the instrument through an arc, because if the 
instrument is moved too rapidly the maximum telemetry data 
rate is exceeded and erroneous angle data will be trans- 
mitted. This can only be detected by checking the original 
initial peinting. Although this problem was observed during 
a manufacturer's demonstration, it did not occur during the 
experimental field work. Pinally, th2 transmitter range of 
5 Km is rather short for the distances used by NOS, which 
can be up to about 10 km. The manufacturer has stated that 
the system range can be easily extanded by increasing the 
transmitter power [Apsey, 1983]. Although this thesis meas- 
ured Aztrac error at a maximum distance of 3 km from the 
shore station, the conclusilosns stated here should not be 
blindly extrapolated to increased ranges. Still, if the 
Aztrac pointing error standard deviation is reduced to its 
angular resolution (0.01 degree) at the very slow angular 
speed of the vessel at long ranges, the error appears +o 
remain acceptable. For example, 0.01 degree pointing error 
at a range of 10 km results in an 2rror of only 1.74 meters. 

The advantages of the Aztrac clearly outweigh its disad- 
vantages, soit is therefore recommended that the Odom 
Aztrac system be inccrporated into the computer based equip- 
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ment used by NOS. [ft should be used only in situations 
where it meets reguired positioning standards based on the 
value for gd, found in this investigation. It is also recom- 
merded that the systematic error induced by the time lag 
discussed here should be accounted for in operational stan- 
dards for range-azimuth positioning using the Wild T-2. 
This could be done fer semi-automat2d systems by providing 
an automatic cadio signal to the observer on shore that 
precedes the actual depth measurement by one second. A 
Simpler method of reducing this systematic error could be a 
limit on the speed of the vessal, depending on the accuracy 
standard required for the particular survey. 

A final recommendation must be nade regarding the posi- 
tioning standards of NOS. At pres2ant there are conflicting 
standards for a given survey scale, depending on whether 
electronic, hybrid (inciuding ranga-azimuth), or visual 
methods are used. For example, positions for a 1:5,900 
scale survey may be required to have an accuracy of either 
Heoem@ecers d.~., OF 2.5 %82ters at some unspecified prob- 
ability, depending on whether microwave ranqe-range or 
range-azimuth methods are employed. The NOS is certainly 
the most progressive hydrographic organization with regard 
to poSiticn error specifications, but it is recommended “hat 
Pmemconcept of probability be applied to all positioning 
methods and not only electronic onss. FPurther, if the 
meaning of "saldom exceed" in the THO standards is to be 
interpreted as a 90% probability circle, then all the NOS 


standards should be changed to steflact that standard. 
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GEODETIC 


Station Names 


SOPAR (1947) 
USE MON (1978) 
MUSSEL (1932) 
AZTRAC 

T2 

GEOCEIVER 
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APPENDIX C 


EMPIRICAL PROBABILITY DENSITY FUNCTION PLOTS 
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Figure C.1 Probability Density Plot: 300 m arc. 
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PaeguceuC., 2 Probability Density Plot: 500 m arc. 
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Pigewee Cc, 3 Probability Density Plot: 700 m arc. 
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Figure C.4 Probability Density Plot: 1000 m arc. 
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Figure C.5 Probability Density Plot: 1500 m arc. 
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Figure C.6 Probability Density Plot: 3000 m arc. 
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